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ABSTRACT ey ee nerd 

This thesis presents a preliminary statistical analysis 
of the daily and weekly rainfall for the Monterey Peninsula, 
in central California. The analysis begins by examining 
the daily rainfall data, also the relationship among the 
length of the storms, amount of rainfall in the storms and 
length of the successive days of rain. Also included is a 
study of the distribution of the amount of rainfall in the 
storms. Also study of the distribution was carried out for 
non-zero weekly rainfalls. 4x4 contingency tables are used 
to identify dependence/independence among the weeks ina 
given month. Also, 2x2 contingency tables are used to 
examine dependencies between weekly rainfalls; logistic 
analysis is used as a parametric model for dependence. 

This paper attempts to analyze rainfall data in the 
Statistical sense. No attempt is made to provide a physical 
explanation of the findings from the point of view of a 


Meceorologist. 
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fee NE RODUC TION 


The Monterey Peninsula Water District, in the central 
California coastal area, has as one of its responsibilities 
the duty to recommend and/or impose water rationing on its 
constituents. 

To do this in a rational way requires the district to 
have some way for predicting future water availability. This 
thesis analyzes rainfall data for the Forest Lake station of 
Monterey by purely statistical methodology in order to iden- 
tify possible ways for predicting future water availability. 

No strong evidence for useful procedures has been uncov- 
ered in this thesis, although some week indications of possi- 


ble dependencies have been found. 
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It. DATA 


A. GENERAL 

Daily and weekly rainfall data were used in this analy- 
sis. The data were accumulated at the Forest Lake station 
of Monterey, in Central California. Rainfall data has been 
gathered by the California American-Water Company since 1891. 

Although this data set started guite early, the data 
prior to 1938 has frequent missing observations. Therefore 
this data set includes observations from January 1938 through 
December 1974. Appendix A contains a listing of the daily 
memtall data. 

Weekly rainfall amounts have been obtained by summing 
daily rainfall amounts, starting from the beginning of Octo- 
ber and running to the end of the April (which is considered 
the rainy season for Monterey Peninsula area) for the 36-year 
period (1938-1939 through 1973-1974). A week is defined in 
terms of Julian dates rather than the usual calendar week. 
For example the first week of January is defined to include 
the first seven days of the year. For other definitions of 


weeks, see the Table of Symbols. 


B. WEEKLY DATA 

Appendix B contains a listing of weekly rainfall data. 
Appendix C shows plots of the weekly rainfall. As can be 
seen the data are strongly seasonal. This is enough to 


indicate that it is quite non-stationary. 
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Means and variances of the weekly rainfall are shown in 
Tables 1 and 2 for weeks with and without positive rainfall 
respectively. Figures 1 and 2 show plots of means and vari- 
ances for weeks with and without positive rainfalls. On 
the figures week 1 represents the first week of October as 
Ol and week 31 represents the end of April and beginning of 


May as AM. 
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III. ANALYSIS OF STORMS 


A. GENERAL 

This analysis 1S carried out on rainfall data from Octo- 
ber to May for a 36-year period (1938-1939 through 1973-1974). 
The rainy period is considered to run from October to May 


for the Monterey Peninsula. 


B. STORMS 

It will be said that a storm lasts exactly n days if 
there are exactly n consecutive days having rainfall greater 
than 0.02 inches. For example, if on January 10th there is 
no rainfall and on the llth, 12th, and 13th of January there 
@eenwo.s0", 0.15", 1.15" recorded in inches of rain respec~- 
tively, and on the 14th of January there is again no rainfall, 
this means that a storm of length of duration 3 days has 
eccurred, and the amount of rainfall in this storm is 1.60" 
meso + 0.15 + 1.15 = 1.60). 

Based on the above definition the historical lengths of 
Storms and the amount of rainfall in the storms will be 
examined. 

Appendix E shows the histograms of length of the storms, 
as denoted by LS, amount of the rainfall in the storms, as 
denoted by AR, and length of the non-rainy period after the 
storms, as denoted by LN for October through April and De- 
cember through February in the 36-year period. The rainy 


period in the December through February months is more 
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memegencous than the October through April period so for this 
reason the December through February period is also examined. 
Figures 3 through 5 show the time series plots of the LS, AR, 
and the LS in the 36-year period. On the figures, dot (-) 
moaacates the beginning of each year from 1938-1939 through 
1973-1974. Figures 6 through ll show the plot of LS against 
AR, LS against LN, and LN against AR for the October through 
April and December through February months in the 36-year 
period. 

Plots of LS versus LN indicates that if the length of 
the storms is small then the following non-rainy period 
(dry period) is large. Table 3 shows the LS and the mean 
of the length of the following non-rainy period for the Octo- 
ber through April and December through February months. 
Figure 12 shows the LS versus the mean of the following non- 
rainy period for the October through April and December 
through February in the 36-year period. The length of the 
storms (LS) and amount of the rainfall (AR) in the storms 
appear to have a linear relationship. By using median methods 
the slope of the length of the storms against the amount of 
the rainfall in the storm was computed as 0.71 for October 
through April and 0.79 for December through February. 

Here are some statistics from the historgrams of LS, 
AR, and LN (Table 4). 

Appendix E shows the histograms of the amount of rain- 


falls in exactly n days lasting storms which are made for 
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Table 3: MEAN CFP LENGTH OF THE NON-RAINY PERIOD FOLLOWING 
EXACTLY N DAYS LASTING STORMS IN THE 36-YEAR PERIOD. 


OCT-APR DEC-FEB 
STORM MEAN OF LN MEAN OF LN 

IN DAYS IN DAYS IN DAYS 
LS=1 TEBE 6.04 
LS=2 6.71 5.92 
LS=3 6.47 5.34 
LS=4 6.66 3.83 
LS=5 6.44 6.17 
LS>6 5.91 3.76 
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Figure 12. 
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Plot of the LS against mean of the non- 
rainy period for October-April and 
December-February 
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maple 4: SOME STATISTICS FOR THE LS,AR, AND LN 


OCTOBER-APRIL DECEMNBER-FEBROARY 

Min. Mean Max. Min. Mean Max. 
LS in days We 2.09 14, 1; 2.37 — 18: 
AR in inches 0.03 0.71 7.86 0.03 0.83 7.86 
LN in days tie 6206 57. Ue 504) 745. 


feel, 2, 3, 4 days in the 36-year period for October through 
April and December through February. 
From the total number of occurrences for each exactly 
n days lasting storms estimated of occurrences for a year 
were computed in the following way and are shown in Table 5. 
This suggests that many of the small storms, especially 
exactly one day lasting storms, occur during the months of 
October, November, March, and April and many ofthe large 


storms occur during the months of December, January, February. 


Sg, 





ieee 3:5 ESTIMATED + OF OCCURRENCES/YEAR FOR THE EXACTLY 
N DAYS LASTING STGRMS IN THE 36-YEAR PERIOD. 


OCT ,NOV,MAR, APR DEC,JAN, FEB 


STORMS ESTIMATED ESTIMATED 


(EXACTLY) OCC URRENCE/YEAR OCCURRENCE/YEAR 


Ls=1 6.52 4.64 
LS=2 2.97 So 
LS=3 1. 06 Weck 
LS=4 0.78 0.84 
LS=5 0.28 0.47 
LS>6 0.36 0.53 
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IV. ANALYSIS OF WEEKLY RAINFALL 


A. DISTRIBUTION OF POSITIVE WEEKLY RAINFALL 
me Lheory 

In this section we will explore the distribution of 
positive weekly amounts of rainfall. A week is said to have 
a positive amount of rainfall if the amount of rain exceeds 
0.02 inches. Table 6 1s a listing of estimates of means 
and variances for positive weekly rainfalls. They were com- 
puted by using the program "HISTG". Figure 13 shows a plot 
of the weekly means and variances respectively. The pattern 
of weekly means again indicates nonstationarity of rainfall 
with more rain on the average falling during the weeks D4, 
J2, J3, JF, F3, FM, M4 where the average amount of rainfall 
in these weeks is more than 1.00 inches. 

eset oical Distribution of Positive Weekly Rainfall 

Pieoncer tO fe a theoretical distribution to the 
positive weekly rainfall it is necessary to use some sta- 
tistics computed from the data such as the mean, median, 
variance, standard deviation, quartiles, coefficients of 
variation, and skewness and kurtosis coefficients. Agreement 
of these statistics with the theoretical values for the ex- 
ponential model (or some others) helps to identify and sup- 
port that model. 

Table 7 shows a listing of means, standard deviations, 


and coefficients of variation of the non-zero (positive) 


Ce, 








MEANS AND VARIANCES FOR POSITIVE WEEKLY RAINFALL 
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Figure 13. Weekly means and variances for the positive 
weekly rainfall data in the 36-year period 
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AND COEFFICIENTS 


STANDARD DEVIATIONS, 
OF VARIATION FOR POSITIVE WEEKLY RAINFALL. 


MEANS, 
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weekly rainfalls. Figure 14 shows the weekly coefficients 
of variation for the positive weekly rainfalls. 

As can be seen, the means and standard deviations of 
the data are approximately equal to each other, and the 
coefficients of variation of the data are close to one. 

These facts indicate that the distribution of the data may 

be approximately exponential. This is suggested since for 
data from an exponential distribution, the mean of the data 
equal the standard deviation of the data, and the coefficient 
of variation is one. 

The fit of an exponential model for the positive 
weekly data will not be explored in more detail. 

A plot of the means against standard deviations of 
the positive weekly rainfall is shown in Figure 15. These 
data appear to have a linear relationship. By using the 
median method (McNeil [Ref. 3]) the slope of the means versus 
standard deviations was computed as 0.84. This indicates that 
the means of the data are little higher than the standard 
deviations of the data. 

Next, by using the subroutine EXPLT (which is a 
NONIMSL subroutine at the Naval Postgraduate School computer 
library), a plot exponential scores versus observed scores of 
the data can be examined. If the data is distributed exponen- 
tially, the the plot should be a straight line. Appendix E 
shows the exponential scores versus observed scores for the 
positive weekly rainfall. Most of the weeks seem to fit the 


exponential. 
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COEFFICIENT OF VARIATION 


ae oe 1G) ea Awe 2515 SB) Sieye 


WEEKS OF THE DATA 


Figure 14. Weekly coefficients Of variation for the 
positive weekly rainfalls in the 36- 


year period 


44 





potzed zeak-9¢ 


S94} UT [TTejyutTer ATYSOM SaTATSOd 3uyA ACF 
SUOTReTASp prepueAsS yYsutebe sueow JO ROTd “GT eanbty 


eT°O 


6€°O 


6$°0 


08°0 


00°T 


O2°T 


lon 


SLINA T0-3602°O =new FFIIVIS-A 


SUCOW SLINA TO-3T9TPO enon 8 831V9S-K 
6g °T €z°T 06°90 8¢°0O 92°0 
SHFSHSSHSSHEFESARSESSRSHCSSESSSSSES SSS TS SHS SSEHSSSCS+SSaseSs ees? Sess SHSTH+ OS SSSS SHS + SOHO SSESSES CE 
+ ry 
s ® + 
¢ s 
* e 
* e 
td ¢ 
s @ 
* * * 
e as + 
s ; i @ 
+ < ry 
” ? e 
* = 5 e 
* e 
* ss ¢ 
bd e 
* e 
td ¢ 
* ie 
* * 
4 ® ® + 
e x4 ~ 
* . « 
s = es 
” * es 4 
s ¢ 
’ ¢ 
] a “t 
¢ ¢ 
* a 
* * 
* * 
* sd ; * 
e * 
* * 
¢ e 
* ~ Ps 
s f e 
s # 
. * 
+ td + 
* s 
s ¢ 
s as e 
* m 2 e 
e . e 
t s 
s e 
* i = ¢ 
* td 
+ > 
s td 
¢ e 
s * 
* “d * 
e e 
s * 
* é 
e * * 
. e 
* °* > 
*eett* SHOTS SH OST HEE SHE HE+EHSEHESCHEH+E HS SSEES aera en  ee eee 776 
SS° €z2°l 6° S° 92° 


etT°O 


6€°O 


6$°9 


03°0 


00°T 


Oz°T 


wet 


TAS 
“p3As 


45 





Also, estimated and theoretical values of median, 
lower and upper quartiles were examined against each other. 
Table 8 shows the estimated and theoretical values of the 
median, and lower and upper quartiles for the positive weekly 
rainfalls. Figures 16 through 18 show the plot of estimated 
values against theoretical values. They seem to support the 
exponential model for the positive weekly rainfall. Again 
using the median method the slope of estimated values against 
theoretical values were computed and are shown in Table 9. 

A 90% confidence interval for the average positive 
weekly rainfall can be computed by assuming exponential dis- 
tribution for the weeks which are shown in Table 10. Let 
X be average positive weekly rainfall with n observations, 
then the formula for (l-a)% confidence interval for the mean 
with a = 0.10 is 


2 
feet t=—a/2) | 


Zz 
ee 2) 
ora Ss) Bice = (nla 


Zi 


aT 


B. SKEWNESS AND KURTOSIS ANALYSIS FOR POSITIVE WEEKLY 
RAINFALL 


1. Theory 
In general the third moment of a distribution is con- 
Sidered to be a measure of skewness. If the distribution of 
a sample is symmetric, its third moment abaut the mean will 
be zero. If the distribution is skewed to the right, the 
third moment about the mean will have a positive value, be- 


Cause the large size of observations on the long tail will 
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LOWER 


ESTIMATED AND THEORETICAL VALUES OF MEDIAN, 
AND UPPER QUARTILE FOR THE POSITIVE WEEKLY RAINFALL 


Table 8 


THEOR. ESTIM. THEOR. 


ESTIM. 


LOWER LOWER UPPER UPPER 


THEOR. 
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WEEK MEAN MEDIAN MEDIAN QUARTILE QUARTILE QUARTILE QUARTILE 
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Table 9: 


Bol « 


EST. 


EST. 


SLOPES OF THE SSTIMATED VALUES AGAINST THEORETICAL 


VALUES FOR MEDIAN AND LOWER AND UPPER QUARTILES. 


SLOPE 
MEDIAN YS THEO. MEDIAN 0.96 
L.QUAR. VS THEO. L.QUAR. 0235 
U.QUAR. VS THEO. U.QUAR. 0.88 


3 dt 


ENDER CERES 





90% CONFIDENCE INTERVAL FOR AVERAGE POSTIVE WEEKLY 


RAINFALL WITH EXPONENTIAL DISTRIBUTION 


Table 10 


GPPER LEM 


LOWER LIAIT 


STD DEV 


MEAN 


WEEK 
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more than offset the greater number of smaller observations 
on the shorter tail of the distribution. Hence for a posi- 
tively skewed distribution (1.e., one with the long tail to 


the right), the third moment (u.,) will be positive. For 


23 
these reasons the third moment around the mean is taken as 
a measure of the absolute skewness of a distribution. The 


theoretical skewness for a random variable X is given by 


, eg EL(XHE(X))°] 
S (var [x])?/? 
The empirical skewness for data {X., 1 = eee 
having average x is given by 
n 
= ) (x. -x)? 
ee! 
1 n 
it —2,3/2 


If a distribution 1S symmetric, will be zero. If 


er 
the distribution of X 1S exponential with parameter A = l, 
then Y1 has the value 2 (see Appendix D for the algebraic 


Gomputation). 


The kurtosis for a random variable X is given by 


BL(X-E(X)) 4) _ 


: 
Z [var (Xx) ]? 


The empirical kurtosis for data {X51 ite lene Deh mang 


average x is given by 


Sye 





= | d, (X, -X) =) 
= 


n 
a Y (x, -x) 1° 
i=l 


If the fourth moment (uy) about the mean of the 
random variable X is large relative to the variance cane 
it indicates relatively large tails. 

For the normal distribution, u 4/0" has the value 3. 
Since the normal distribution arises very frequently and is 
often used as a basis of reference for distributions that are 
not normal, the quantity Y9 is defined so that it will be 0 
when a distribution has the kurtosis of a normal distribu- 
eeen. Thus og O means that a distribution has a sharper 
peak, thinner shoulder, and fatter tails than the normal dis- 


mialoution. < 0 means that a distribution has a flatter 


2 
peak, fatter shoulders, and thinner tails than the normal 
Seeeettoution. Cramer [Ref. 5] and Duncan [{Ref. 6] contain a 
discussion of the skewness and kurtosis coefficients. 

If the distribution of X is exponential, then Y5 
has the value 6 (see Appendix D for algebraic computations). 
Appendix D also contains a discussion of the sample proper- 
ties of the skewness and kurtosis coefficients. 

It 1s suspected that the sample size has some effect 
on the values of sample skewness and sample Kurtosis. To 


study this effect a simulation study was done as described 


below. 
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By using the random number generator (LLRANDOM) N 
independent unit exponential random numbers were generated 
as a sample. Then the sample skewness and kurtoSis were 
computed from the sample. M independent replications of the 
procedure were done and sample means and standard deviations 
were computed for the skewness and kurtosis. 

Appendix D shows the simulation results for sample 
skewness and sample kurtosis of a unit exponential distribu- 
tion with various sample size (standard deviations of them 
are given in parentheses). 

As can be seen in Appendix D if the sample size 1s 
small, then the sample skewness and kurtosis values are 
smaller than their theoretical values. When sample size is 
between 2000-3000 they reach the theoretical values 2 and 6 


for the skewness and kurtosis respectively. 


2. Skewness and Kurtosis Analysis fon Positive Weekly 
Rainfall 


Histograms, and plots of exponential scores versus 
observed values of positive weekly rainfalls are given in 
Appendix F; they suggest that the distribution of weekly 
positive rainfall is approximately exponential. The examina- 
tion of classical skewness and kurtosis coefficients of posi- 
tive weekly rainfall can be used to further examine the fit 
of this exponential model. 

Table 11 shows a listing of the values of skewness 


and kurtosis for the positive weekly rainfall. 
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SKEWNESS AND KURTOSIS FOR WEEKLY RAINFALL DATA 


Table 11 


NUMBER OF YEARS 


POSITIVE RAINFALL 


NORE oo to 


SKEWNESS 


WEEK 
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Figure 19 shows the skewness and kurtosis of the 
ositive weekly rainfall, and a plot of skewness against 
urtosis values is shown in Figure 20. 

Figure 19 shows that the skewness and kurtosis values 
re reasonably stable. But relatively high values occur early 
n the year. There is some tendency for the distribution to 
shange throughout the rainy season. Average skewness and 

curtosis values of the positive weekly rainfall are shown 


in Table 12. 


Table 12: ESTIMATED AND SIMULATED MEANS POR THE 
SKEWNESS AND KURTOSIS COEFFICIENTS 


SKEWNESS KURTOSIS 
ESTIMATED MEAN e239 1.19 


SIMULATED MEAN 1.45 1.34 


The estimated kurtosis is somewhat lower than its simulated 
yalue. This might suggest that the distribution Ot POSsLeLve 
weekly rainfall has a shorter right tail than the exponential. 
Perhaps a Gamma or Weibull distribution would fit these data 
slightly better. 

Table 14 shows the estimated and simulated values of 
the skewness and kurtosis coefficients for the actual sample 
sizes observation. By using the median method the slopes of 
estimated values against simulated values were computed for 


skewness and kurtosis and are shown in Table 13. 


c7 





Neary ShenNess 


oil ell allo a bass 


WEEKLY 





Jaa Coeen es iimc DATA 


= SS ———— oo 


Ie 


| 


: | 
0. : ji MA WA AL \y 


life Oe lhe BATA 


Figure 19. Weekly skewness and kurtosis coefficients for 


the positive weekly rainfall 
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EST. 


ES. « 


Table 13s 


SLCPES FOR SKEWNESS AND KURTOSIS 


SLOPE INTERCEPT 


SKEWNESS VS SIM. SKEWNESS 0.19 1.20 


KORTOSIS VS SIM. KURTOSIS 0.40 1.32 
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ESTIMATED AND SIMULATED VALUES FOR SKEWNESS 


Table 14 


KURTOSIS FOR SAME SAMPLE SIZES. 
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As can be seen they are not fitted very well. Esti- 
mated values are lower than the simulated values for the 
skewness and kurtosis. Figures 21 and 22 show the estimated 


values against simulated values for the skewness and kurtosis. 


fee 2X2 CONTINGENCY TABLES 
1. Summary 

The idea to be explored in this section is whether 
or not some weeks of the weekly rainfall, to be called the 
control, may be used to predict in some way the behavior of 
another weeks of the weekly rainfall, to be called the 
complement. 

Let X be a week of a month, to be called the control, 
and let Y be a week of a month, which is different from X, 
to be called the complement. It is necessary that Xn Y = Q; 
that is, the intersection of these two weeks is empty. The 
weeks are then compared for some quality in X and for some 
quality in Y. The question is then; does the presence (or 
absence) of quality in X affect the presence (or absence) of 
quality in Y? An example of a typical table is shown below 
maeraigure 23. 

The elements, a typical one of which is Ba 6 repre- 
sent the number of years which display quality 1 in the con- 
trol and quality j in the complement. The marginal entries 
ns and n 4 represent the number of years for which the con- 
trol has quality i and the number of years the complement 
has quality j respectively. The overall number of years, N, 


appear in the lower right of the table. 
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megure 23. Typical 2X2 contingency table 


Conover [Ref. 2] gives a discussion of the theory and 


use of contingency tables. Let re be the probability that 


any given year will have a control quality i and a complement 


quality j. Then estimates of the ogee are 
ye = ad 
a3. = nN 
o = Ba 


If the control and complement are independent, 


835 = La a pei “ald 3. 
These simple assumptions allow for an investigation 
of the possible interrelationships between the weekly 
rainfalls. 
The question of independence may be approached in 


the was as described below. 
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mmeecnl-oquared Test for Independence 


Let Bia equal n, n ad he ici ea Orucm eee CoML mgency 


table the test statistics are given by 


2 2 (n, .-E,.)° 
a ee 
tel 4e0 ij 


Seesamplifying for the calculation 


2 Z n.* 
Q= ) J -l-n 
i=l j=l “ij 


Mae exact distribution of QO is difficult to tabulate 
because of all the different combinations of values possible 
mor Mii. Myo7 Aoy1 Agd- Therefore the large sample approxima- 
memeis used for the distribution of Q; this turns out to be 


chi-square distribution with one degree of freedom. 


Hypothesis 

Hy? O55 = Jie ; FOR ae. 2g 

Hy: i — - ean - for some 1,}) 
ala) 5 ee as 


Decision Rule: Reject Ho is Q exceeds the (l-a) 
gQuantile of a chi-square random variable with one degree of 
Freedom. The approximate level of significance is then a. 

Figures 24 through 28 show 2x2 contingency tables for 
the weekly rainfall. In these tables, X = 1 occurs when the 


week i falls below its mean of weekly rainfall and X = 2 
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28 8 36 
g=0.32 
(a) 
03 VS O88 
1 2 
—  «t 
1 18 7 25 1 
} 
. | Z j 11 2 
EE ee | 
27 9 36 
Q=0.39 
{C) 
ON ¥S N1 
1 2 
( 
Z a ! 2 2 2 
24 V2 36 
Q=0.67 
(e) 


Figure 24: 2x2 contingency tables for 
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O72 VSe03 


25 14 36 
Q=1.83 
(D) 
O4 VS ON 
7 2 
| 21 | 6 |27 
i. 2 4 | 9 
26 10 36 
Q=1.66 
(d) 
N1 VS N2 
1 2 
2 
| bs) y | 12 
(272 14 36 
Q=2. 86 


(2) 


weeks 01,02,03,04,0N,N1. 





Me ¥S N3 N3 VS NG 











1 2 1 2 
1 | 19 : 3 | p>? 1 22 6 | 28 
2 9 < 14 2 5 | 3 8 
ee SS eee 
28 8 36 a7 9 36 
Q=2.41 0=0. 86 
(a) (b) 
NG WS D1 D1 VS D2 
1 2 
1 16 11 aT 1 23 
a te 
2{ 7 | 2 | s : 13 
23 13 36 25 11 36 
Q=1.00 QO=0.54 
(c) (d) 
D2 YS D3 D3 VS D4 
1 2 1 2 
1 16 9 25 1 16 4 20 
: —s ae eee ees 
p | 4 | 7 11 2 {| 8 | 8 t 16 
i Se a ee 
20 16 36 2u 12 36 
Q=2. 36 0=3.60 
(2) (f) 


Pigure 25; 2x2 contingency tables for weeks N2,N3,N4,D1,D2,D3,D4. 
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D4 VS J1 oa VS 708.2 
1 Z 1 2 


ae 


24 1 16 7 Z3 


— | | 
14; 18 { 6 { 
2 | 5 | qT 2 2 7 6 | 13 
CO re 
23 13 





36 23 13 36 
Q=3.85 Q=0.89 
(a) (b) 
J2 vs 33 J3 vs Ja 
1 2 
1 | 18 mn } 22 
a 2 4 
Pa) Sohal cee 
24 12 36 
QO=5. 84 
{d) 
JF VS F1 
1 2 1 2 
—_ ff °°  . ri | [ae 
1 14 10 24 1 16 7 | 23 
+4 
2 | g | 3 12 2 | 9 Z 4 | 13 
23 13 36 25 11 36 
Q=0.96 © 9=0.01 
(e) (f) 


Figure 26: 2x2 Contingency tables for weeks D4,J1,J2,J3,J4,J0F,F1. 
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mys F2 R2 VS F3 





1 D 1 2 
1 | 17 | 8 | 25 1 16 6 22 
[pf aa. 
2 5 6 | (Pa 2 7 7 | 14 
es 
20 14 36 23 13 36 
0=1.63 Q=1.92 
(2) (b) 
F3 VS FM PM VS M1 
1 2 1 2 
1 17 6 | 23 1 i 17 | 8 25 
2 | 8 | 5 13 2 | 5 ) 6 11 
ee | eee | 
25 11 36 22 14 36 
9=0.60 9=1.63 
(<c) (d) 


M7 VS M2 2S a8 


| | | 
Z | 2) | 5 14 Z | 6 7 | 13 


23 13 36 23 is 36 


Pigure 27: 2x2 Contingency tables for weeks F1,F2,F3,FN,41,42,M3. 
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M3 VS M4 MG ¥S At 


1 2 1 2 
—__ a eT ae 
1 16 | 7 23 1 22 5 | 27 
| 
1 
2 1 2 | 13 2 f 4 5 | 9 
27 9 36 26 10 36 
Q=1.00 Q=4.62 
(a) (b) 
Al VS A2 A2 VS A3 
1 2 1 2 
— t 
1 20 6 26 1 21 | 4 | 25 
| ) 
2 | 5 | 10 2 7 | 4 | 11 
0 eo ee ee 
25 11 36 yr: 8 36 
Q=2.47 Q=1.83 
(Cc) (d) 
A3 VS AY AY ¥S AM 
1 2 1 2 
1 | 17 | 11 28 1 | 16 ! 7 | 23 
| | 
m cs (| | 2} oe | 5 1143 
www a | ee 
23 13 36 24 12 36 
Q=0.55 9=0.24 : 


Figure 28: 2x2 Contingency tables for weeks M3,M4,A1,A2,A3,A4, AM. 
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occurs when week 1 falls above its mean of weekly rainfall. 
Similarly, Y = 1 occurs when the week itl falls below its 
mean of weekly rainfall and Y = 2 occurs when the week i+l 
falls above its mean of weekly rainfall. 

The control consists of the binary category of rain- 
fall for week i and the complement consists of the binary 
category of rainfall for week i+l in the 36-year period. 

Results of the test statistics suggest that the weeks 
Nl vs N2, D4 vs Jl, J3 vs J4, M2 vs M3, and M4 vs Al are not 
independent with a = 0.10 significance level. The rest of 
the weeks appeared independent of each other at a = 0.10 
Significance level. Table 15 shows the test statistics and 
their chi-square values with a Significance levels. 

The expected number of 30 Ve random variables that 
would be over .10 is 3. So it seems that the result could be 
explained by random sampling. 

Figures 29 through 33 show the 2x2 contingency tables 
for which the weekly X and Y values are determined differently 
from those shown above. On these figures, X = 1 occurs when 


the week i has no rainfall (no rainfall means rainfall of 


less than 0.02") in year t and X = 2 occurs when week i has 
Painfall in year t. Similarly, Y = 1 occurs when the week 
i+l has no rainfall in year t and Y = 2 occurs when the week 


i+] has rainfall in year t. 
The control consists of the binary category of rain- 
fall for week i and the complement consists of the binary 


category of rainfall for week itl in the 36-year period. 
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Table 15: TEST STATISTICS, APPROXIMATED CHI-SQUARED 
AND SIGNIFICANCE LEVEL FOR THE YEEKS. 


TEST APPROXIMATED § SIGNIPICAiCE 
WEEK § STATISTICS x4) VALUE | LEVEL 
01 VS 02 0. 32 ae 0.60 
02 VS 03 1.83 —e 0.18 
03 VS o4 0.39 oe 0.53 
04 VS ON 1. 66 ay 0.20 
ON VS NI 0.67 ke 0.43 
N1 VS N2 2.85 an 0.09 
N2 VS N3 2.41 ‘oe ene 
43 V5 NU 0. 86 ~~ 0.37 
N4 VS D1 1.00 -— 0.32 
D1 VS D2 0.54 ee 0.48 
D2 vs D3 De BAe X99 | 0.13 
D3 VS DY 3. 60 ey, 0.06 
DU VS J1 3. 85 ee 0.05 
mys J2 9.89 ner 0.35 
J2 vs 33 0.45 ee 0.51 
33 YS Ju 5. 84 > 0.02 
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om VS O02 
1 2 
—___ 

1 | 13 | 14 | BY 1 
2 | 6 | 3 t 9 2 
———e 

19 17 36 
Q=0.93 
(a) 
03 VS 0&8 
1 2 
1 12 8 20 1 
mes) || 7 jit 2 
es | 
21 Ls) 36 
Q=0.05 
(Cc) 
ON VS Nl 
1 2 
1 11 11 | 22 1 
2 5 3 | 14 2 
Eee 
16 20 36 
Q=0.71 
(e) 
Figure 29: 


ces 


O20 75: 103 
1 2 


11 8 19 
| | 
Pe), 
20 16 36 
Q=0.09 
(b) 
04 ¥S ON 
15 6 21 
17 | 38 its 
22 14 36 
Q=2. 26 
(d) 
NI VS N2 
1 2 
5 11 16 
re a 
3. | 17 | 20 
a | 
8 28 36 
Q=1. 36 


(£) 


4 


2x2 contingency tables for weeks 01,92,03,04,0N,N1. 





N2 VS N3 N3 VS NY 





1 2 1 2 
1 | 5 | a | 8 1 | 8 : 9 | 17 
21 12 #+| «6 | 28 Dn ener eae fi) 19 
Oe EE es | SS ee 
17 19 36 14 22 36 
Q=0.96 Q=0.91 
(a) (Db) 
NG VS D1 D1 VS D2 
1 2 
1 14 1 | 3 5 | 8 
2 22 2 11 17 28 
36 rT ey 
Q=0.54 Q=0.01 
(Cc) (d) 
D2 VS D3 D3 VS D4 
1 2 1 2 
——_ ., 9 Pao -— ~ a 
1 4 10 14 1 | 5 | 3 | 8 
2 4 18 | 22 2 | 7 21 | 28 
28 36 12 24 36 
Q=0.54 Q=3.98 


(e) (f) 


Pigure 30: 2x2 contingency tables for weeks N2,N3,N4,D1,D2,D3, D4. 
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DY VS Jl 


1 2 
1 | 3 | 9 | 12 
p 7 17 | 24 
10 26 36 
Q=0.07 
(a) 
mo ¥S 33 
1 2 
1 : 8 | 3 11 
| 
mr 7 «6| «(1s (| (2s 
i 
15 21 36 
Q=6.29 
(Cc) 
Ja vs JP 
1 2 
1 3 6 9 
2 4 23 | 27 
29 36 
Q=1.483 


Figure 3l: 2Zx2 Contingency tables for weeks D4,J3J1,J2,J33,J4,J0F,F1. 


Oi Misa 2 


Jo 21 36 
Q=5.66 
(b) 
J3 VS gu 


| Se 10 | 1S 
4 | 17 f 21 
27 36 
Q=0.95 
(d) 


OF Vs emt 
1 2 


5 2 7 


| u 25 29 
25 36 
Q=9.99 
(f) 





ma vs F2 Pe 7%5 03 





1 2 1 2 
-_— ). <a | 
7 | 2 | 7 i 9 1 yu 5 ; 9 
| | 
7 
2 7 20 {2 2 11 16 {| 27 
9 27 36 15 21 36 
Q9=0.05 Q=0.04 
(a) (b) 
F3 VS PM FM VS M1 
1 2 
| 
1 4 i 11 | 15 6 
2 | 2 ' 19 j 21 30 
aia ca 
6 30 36 36 
Q=1.85 Q=0.23 
(Cc) (da) 
M1 VS M2 M2 VS M3 
1 2 2 
1{o {| a | « zi os |e 
| | [ 
2 | 4 24 | 28 i 20 { 28 
FEE Ee 
4 32 36 25 36 
Q=1.29 Q=0.23 


(e) 


(f) 


Pigure 32: 2x2 Contingency tables for weeks F1,F2,F3,FM,N1,M2,M3. 





M3 VS M& M4 VS Al 


1 2 1 2 
1 | 6 | 5 | +1 1 6 7 | 13 
| | | | 
Z 7 18 1-25 iz ee | ZS 
EE ee | 
13 23 36 15 21 36 
Q=2. 33 Q=0.17 
(a) (b) 
A1 VS A2 A2 VS A3 
1 Z 1 2 
Sf. °° 7 ae re 
1 7 8 5 1 11 3 14 
[ | 
2 7 14 {| 21 Z 11 11 | 22 
a eS [| ee Kee | 
14 22 36 Ze 14 36 
0=0.55 Q9=2.94 
(c) (d) 
A3 VS AY AY VS AM 
1 2 1 Ps 
| | | | 
P q j 10 | 14 2 | 9 13 22 
ee CE 
14 Ze 36 18 18 36 
Q=1.03 Q=1.87 
(e) (£) 


Pigure 33: 2x2 Contingency tables for weeks 43,44,A81,A2,A3,A4,AM. 
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fap5te §6 shows the test statistics and their approxi- 
mated chi-square value with a significance level. 

Results of the test statistics showed that the weeks 
Beever D4, Jl vs J2, J2 vs J3, JF vs Fl, and A2 vs A3 are not 
independent with a = 0.10 significance level. The rest of the 
weeks appeared independent of each other at a = 0.10 signifi- 
cance level. According to this procedure, the overall asso- 


Ciation among the weeks is not strong. 


D. 4X4 CONTINGENCY TABLES 


Ll. General 


Figure 34 shows a typical 4x4 contingency table. The 
table element, Sag represents the number of years which dis- 
Play quality i in the control and quality j in the complement. 


The marginal entries n; and n 5 represent the number of years 


COMPLEMENT 


WEEK. (Y) 
1+] 





mugen sa. Typical 4x4 contingency table 
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meee Le TEST STATISTICS, APPROXIMATED CHI-SQUARED 


AND SIGNIFICANCE LEVEL FOR THE WEEKS. 


TEST APPROXIMATED §§ SIGNIFICANCE 

WEEK © STATISTICS 1 _| VALUE LEVEL 
01 ¥S 02 0.93 ue 0.34 
02 ¥S 03 Q. 09 fare 0.75 
03 Vs o4 0.05 rh 0.79 
04 VS ON 2. 26 «87 0.13 
ON VS NI 0.71 <r 0.41 
N1 VS N2 1. 36 ee 0.25 
N2 VS N3 0. 96 67 0.33 
N3 VS NU 0.91 ee 0.34 
Na VS D1 0. 54 ee 0.48 
D1 VS D2 0.01 fe 0.88 
D2 VS D3 0. 54 oe 0.48 
D3 VS D4 3.94 X95 0.05 
Da vs J1 1 OY ar 0.76 
J1 VS J2 5. 66 <r 0.02 
evs J3 6.29 X99 Oa 
33 VS Ju 0.95 a 0.33 
J4 VS JF 1. 48 ae 0.23 
me YS Fi 9.99 ane 0.01 
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for which control has quality 1 and the numbers of years 
the complement has quality j respectively. The overall nun- 
ber of years, N, is in the lower right of the table. 

Conover [Ref. 2] discusses the theory and use of rxc 
contingency tables. 

Z Classification of Rainfall 

The National Weather Association classifies the daily 
positive rainfall in 3 categories, as light rainfall, moder- 
ate rainfall, and heavy rainfall. Any rainfall amount less 
than 0.02 inches is defined as zero rainfall. The observed 
positive rainfall data is ordered from least to highest for a 
given year. Then, the first 1/3 of the ordered positive rain- 
fall data is called light rainfall, the second 1/3 of the 
ordered positive rainfall data is called moderate rainfall, 
and the remaining 1/3 of the ordered positive rainfall data 
is called heavy rainfall. 

To construct a 4x4 contingency table using the weekly 
rainfall data, a similar procedure was used and is described 
below. 

The weekly positive rainfall data was ordered from 
least to highest for a given month in the 36-year period; then 
3 categories of rainfall were defined as described above. 
Additionally, a weekly rainfall of less than 0.02" inches is 
called zero rainfall and is considered as one category so that 
weekly rainfall data was broken down to 4 categories. Thus 


X = 1 occurs when the week i has zero rainfall in year t, 
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7 occurs when week 1 has ii ale lo 
urs when week i has 


urs when week i has heavy ats oem 


ol. 


Table l7 shows the amount 


moderate Paimicall 2neycas or 


t rainfall in year t, 


aqlij sue, SIGENS OE forwras? 


Of rainfall Limits of 5 


-egories of rainfall for the 7-month rainy period. 
Table 17: AMOUNT OF RAINFALL LIMITS FOR 
THE WEEKLY RAINFALL IN INCHES 
LIGHT MODERATE GEAVY 
MONTH RAINFALL RAINFALL RAINFALL 
OCTOBER 9.03-0.12 0.13-0.34% 0.35-0ver 
NOVEMBER 9.03-0.31 0.32-1.00 4,.Q01-over 
DECEMBER 9.03-0.37 9.38-1.07 4.08-over 
JANUARY 9.03-0. 38 0.39-1.34 4.35-0Vver 
FEBRUARY 9.03-0. 40 9.41-1.19 4.20-over 
MARCH 0.03-0. 26 0.27-0.86 0.87-over 
APRIL 9.03-0. 17 0.18-0.690 0.61-over 


ae Chi-Squared Test for Independence 
A 4x4 contingenc 


or rainfall in week i; and the cate 


VacapLe Was construc 


Gory Bye sea geet let rg 


and X = 


tea the category 





9y grouping all weeks in a month together. Table 18 shows 

-he name of the weeks which were included to construct the 
sontingency table for each month. Each month was analyzed 
Sa@ebately Since the rainfall data is seasonal. Each paired 
week contains 36 observations so that the total 144 (4x36 = i144) 
observations are formed in each contingency table. For exam- 
buen the December contingency table the paired weeks used 
g@memtD1,D2), (D2,D3), (D3,D4), (D4,J1). The chi-square test 
was used to explore possible relationships between rainfall 


in successive weeks for a given month. 


Table 18: WEEK NAMES FOR CONSTRUCTING 4x4 CONTINGENCY 


TABLES 
MONTH WEEKS 
OCTOBER 01,02,03,04,0N 
NOVEMBER N1,N2,N3,N4,D1 
DECEMBER D1,D2,D3,D4,J1 
JANOARY J1,3J2,J3J3,J34,J0F 
FEBRUARY Ui eg ec, OO pr 
MARCH PM,M1,82,N3,44 
APRIL A1,A2,A3,A4,AM 


Let Bis equal n, n af Ake Then for a 4x4 contingency 


table the test statistic is given by: 
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ne A (n,.-E,.)° 
QO _ ) ) —=4—=1— 
i=l j=1 ij 


Se simplifying for the calculation 


4 4 aes 
eee eee 
neal Gib ss 


ie exact GlStribution of Q is difficult to tabulate 
because of all the combinations of values possible for ni 3: 
Therefore the large sample approximateion, the chi-square 
distribution with 9 degrees of freedom is used for the dis- 
Mereution of Q. 

Figures 35 through 41 show the 4x4 contingency tables 
for the weekly rainfall. The control consists of the cate- 
gory of rainfall for week 1 and the complement consists of 
category of rainfall for week i+l for a given month, where 
the categories are given in Table 14. 

The results of the test statistics suggest that the 
weeks of January, the weeks of February, and the weeks of 
April are not independent at a= .10 significance level 
Bo, (9) = 14.7 with a = .10). The rest of the months, 
October, November, December, and March appear independent at 
the a = .10 significance level. 

Table 19 shows the test statistics, approximate chi- 


Squared values, and the achieved significance level a. 
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OCTOBER 


as 
; 
os 
|: 
conk 
| 
| ~ 
s 
a 
wa 
co 
—~J 


Nh 
a" 
= 
Lad 
8) 
us 
We) 
UJ 


3 | 10 | > | 4 0 16 
y | 7 | 3 4 4 18 
82 19 21 22 144 


Figure 35: 4x4& contingency table for OCTOBER. 


NOVEMBER 





Figure 36: 4x4 contingency table for NOVEMBER. 
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DECEMBER 


2 a 


| Ce Ts ie ee 
Zz = | 10 [10 10 | 4 | 
3 = at 7 11 | 36 


“ | 17 | | 11 35 
Gy 36 35 144 


Figure 37% 4x4 contingency table for DECEMBER 


JANUARY 


11 | 33 
42 30 35 37 144 


at bat: 22 
of «¢ fe {| ef om | 


Figure 3g 4x4 contingency table for JANUARY 
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FEBRUARY 


3 | 13 | 1 : 7 8 | ao 
A ae ea 
4 6 6 ea 10 : 13 | 35 
a) 37 35 33 144 


Figure 39: 4x4 contingency table for FEBRUARY 


MARCH 
1 Z 3 4 
1 10 | 10 | 2) | 4 29 
| | | 
2 S) ae 2 | 11 | 40 
3{ 13 | 9 | 6 { 10 { 38 
ae ee 


: 4 ) 11 ae 12 | 10 | 37 
6 38 35 


3 35 144 


fe 


Figure 490: 4x4 contingency table for MARCH 
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APRIL 


Oo 
oa) 
On 
a) 


8 | 10 | 31 
68 28 26 22 144 


Figure 41: 4&x4& contingency table for APRIL 
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Table 19: 


PESh SQLALLESL EC, 


APPROXIMATE CHI-SQUARED VALUES, 


SIGNIFICANCE LEVEL FOR THE MONTH. 


MONTHS 


OCTOBER 


NOVEMBER 


DECEMBER 


JANUARY 


FEBRUARY 


MARCH 


APRIL 





TEST APPROXIMATED § SIGNIFICANCE 
STATISTIC CHI-SQUARED LEVEL 
«7-4 (9) , 
Da 1° 49 (9) 0.60 
10.45 «69 (9) 0.32 
8.43 eon) 0.50 
18.64 (799 (9) 0.03 
15.92 X95 (9) 0.07 
11.94 4794 (9) Gaze 
19.71 X99 (9) 0.02 


oe 





SB. LOGISTIC ANALYSIS 
i Lheory 

The logistic analysis to be described in this section 
was developed from Cox [Ref. 1]. 

The basic approach is to view the complement as having 
a binary representation, with success being defined to mean 
that a complement has rainfall, while failure means that the 
Gomplement has no rainfall for week t+l. The problem then 
is to find the conditional probability of a success given 
that the control takes on a particular value. 


Let X, be the rainfall in week t. If the probability 


1 
of “success" = rainfall in week t+l, given Xe 1S written as 
dS Pie aentali in week SPL |G 
the logistic model is 
atBA, 
6 — ae 
tc atBX ; 
l+e 
where the x are the explanatory variables. 
The likelihood function is then 
ate xX es 
EAX,Y3a,B8) = og nttl 
rtiQrer ol a+BX, atBXy 
l+e 1l+e 


where: 


22 





mee weok t+) has rainfall 


t+1 
O if week t+l has no rainfall 


and the log-likelihood is 


N N N atB A, 
Meeere,8) = « ) ¥, +8 } xX ¥, - J In(lte 
t=1 t=L1 el 
The gradient, and Hessians are: 
d9Log L Ps _ — ¢ 
dO ~ a7 te 7 wey ( atEX,) 
- a l+e 
a+B X 
N N © 
dlLog L _ coe <a 
3 7 wy ee L “ak atx, 
Zz ee lt+e 
Hessian: 
ategxX 
N BBs N X.e e 
- ) = 4 ) - X 
ia atgpx atBp 
He = 
atBpxX Op eye 
N X.e E NX. *e ig 
‘ fe 3 ) te 
at ps 2 atBX 
ca (l+e ty t=1 (l+e t,2 
Information Matrix: 
I = -H 
S,S85 IL; 
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to solve for a and g use of Newton's method as follows; 


All necessary elements may be calculated in one pass 
of the computer algorithm. 

One beneficial byproduct of the maximum likelihood 
approach 1s the asymptotic covariance matrix, Fe Sy ae 
Cos [Ref. 3] states that the diagonal elements of this matrix 
provide good estimates of Var(a) and Var (8) under assumptions 
eeenormality. 

Next, by using the diagonal elements of the asymptotic 
covariance matrix to put approximately symmetric confidence 
mmets On a and g; for 90% 

s,s 


Pace eee) le 6 4 (a) 


Pi elG4eNe G63) ae <a + 1.64 VEo"> (8) 


2. Analysis 


In this analysis, given week i has rainfall/no rainfall 
in year t is considered as a control X, then conditional proba- 
bility of success which week itl has rainfall in year t is 


considered as a complement Y. Let 


a if week itl has rainfall 


a 
i+l 
0 if week i+l has no rainfall 
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and 


IL if week i has rainfall 


0 if week 1 has no rainfall 


Then, the model 


atBX. 
PHY Gan, = iL ead 


i 
@ 
i 


atBx, 
1l+e 


and 


g — 
Tag” Pe 


log Odds = ¥ = In( 

Using the computer package "BMDPLR" estimates of a 
and g§ were obtained using maximum likelihood. Other quanti- 
ties computed were conditional observed and predicted. These 
values are shown in Table 20. Also the number of successes, 
number of failures, predicted log odds and parameters a and 
8 and their 90% confidence interval are shown in the same 
tables for the weeks. In the tables weeks are shown with 
their names which are indicated in the glossary of symbols 
tables. 

In the 2x2 tables, the va test results have suggested 
that the weeks 02 vs 03, D3 vs D4, Jl vs J2, J2 vs J3, JF vs 
Fl, and A2 vs A3 are not independent with a = 0.10 signifi- 


cance level. Logistic model also indicates that these weeks 
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foe 20; RESULTS OF THE LOGISTIC ANALYSIS FOR THE WEEKS 


MODEL <: P(02=1{01=x) 


F # OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS O01 
8 7 0% 5333 0.5333 0.134 0 
14 7 0.6667 0.6667 0.693 1 
PARAMETER COEFFICIENT SE COEEJZOS 
Alpha 0.433 0.518 0.7256 
Beta 0.560 0.694 0.806 
Var({ %}=0.15 Var{ 8 jJ=0.65 
90% Confidence limits 
LOwer Aebes 
Beta -0. 762 1.882 
MODEL : 2P(03=1]02=x) 
# OF # OF OBS ERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS 02 
3 11 0.2143 0.2143 12299 0 
11 11 0.5000 0.5000 20 1 
PARAMETER COEFFICIENT SE COREL ZSE 
Alpha -1.299 0.551 -1.995 
Beta 1.299 0.778 1.669 
Var{ . J=0.31 Varf gjJ=0.55 
90% Confidence limits 
LOwer Upper 
Alpha =2en2 2 -0.336 
Beta 0.083 225v9 


3s) 





=X) 


P(O4=1|03 


*e 
= 


HODEL 


\ 
-mM/or 
© 1 

Nw 

Mm) | Or 
OaAtiron 
Mf ee e 
OoOroo 
bf 
Ce 1 
ff 
Hy Pnm 
Oe 1 Wot 
HOiwre 
Ace 1 ne 
MMO « e 
ed ea 
A. 
As 
P| 
et 
Ari 4 
oe 
Pe iwMs 
mratswe 
RIG FUNne™ 
mmm ie »e 
mMmolroo 
Om} 
Os 4 
PrarQ | 
Ontos 
yr 
+ fr, 
Soe | 
hy? | 
OF | 
UO} 
we EAN© 
Sr 
wn | 


PARAMETER 
Ae 
a 
Var{ a] 


Confidence linzits 


Lower 
~0.551 
-Q. 437 


90% 


Beta 


=X) 


P(ON=1]04 


MODEL 


| 
soc 
© | 
| ma 
Cw 1 ow 
OA rom 
“MO * e 
‘ol mele) 
t 
bm | 
AEH | 
Itt 
Het eo 
Coie I ™O 
PM | No 
AcIimMmyN 
MIM 4 « «© 
on coe 
0.04 
A, ft 
bm | 
Fat 
Ate | 
co SS Be Oe ed ©) 
al on I oh 
mOoyne 
hie ye MWYN 
WMI eo 
moO 1oo 
Om | 
Au | 
Peurd | 
OH EDO 
a | 
te Py | 
WN? I 
fx4U? | 
O11 
Ot 
we) 1M 
=) | ie 
Ww 1 


=9.21 


aoeee 
a 
Varf a J 


PARAMETER 
Be 


90% Confidence limits 


er 
s18d 
2-092 


t) 
0 


1{N 1=x) 


252 P NZ 


MODEL 


i 
- | or 
= | 
1own 
wm 1 om 
OA rm 
WHOl ee 
Otlo- 
bm 1 
AE ft 
Ie 1 
Bird § NO 
OH EEO 
HMO i ow 
At | .0cm 
SIM ft © 
mOoOhroo 
0:05 | 
Auf 
4 | 
Ht 
Ari | 
led | iN© 
rH ITO 
fa OQ 4 co” 
PIe<t | WOO 
MMi» « 
mO1oo 
Om | 
Au f 
fer] 
Or Ee 
ect | 
ee ry | 
"| 
fe, 02 | 
Ot] | 
Of 
tO per 
Spicr 
WM) | 


IENT 


OVOor 


furor 
Prue 0 


hoo 


CO 


PARAMETER 
pee 
a 
VarLa J 


90% Confidence limits 


She 
Ey 
PPPS LS, 


J 
7 


Lower 
-90.095 
-0.405 


Seba 


OF 





1{ N2=x) 


P(N3= 


MODEL 


! 
NN} or 
= | 
| =o 
(9K) | ec 
CA IMN 
MNO © « 
O;oo 
] 
bm | 
AR ! 
cH | 
HwI~orsr 
OM 1 Me 
mMmOPere 
Ad I mMW 
EIM 4 2 e 
mohoo 
A:% | 
As | 
bP | 
f 1 
Ak 1 
H41ost 
ie} une 
mmOmeperre 
eT oil Maalte) 
WM e« 
mortroo 
OG } 
O. | 
urd 
OH PEN 
at | a 
1 
mW I 
Ot } 
O } 
tet) FAW 
Se } = 
— | 


ek 
a 
Var{a ]=0.53 


PARAMETER 
Be 


=X) 


T{N3 


—— 
=a 
GP Ob ae ED eee GP <6 ce cae eB ea SS eee ee 


> P(NY 


MODEL 


' 
MIO 
= | 

tom 

ON fem 
QOat rr 
mit = e 
Ohrooe 
an 
AH ! 
BIH} 
Et § aN 
Oe | Oot 
HMI Nw 
Ad mo 
MO ee 
poy adil 
As 
Au } 
Db | 
Fe | 
Ar 1 
dent § oN 
all on OAT 6 
aM Nea 
fd § UNO 
WUMie-e 
mo;oo 
Ot 1 
Qu } 
fort | 
OH 1 Oro 
x | 
te fx; 
wn I 
mn I 
Of] 1 
Ot 
wet 1 nm 
Soi 
wv") 1 


PARAMETER 
ha 


WY) 
rf) 
Ate 
BONO? 
Aare 
riflue ie 
(har Nh mes 
@ 
©) 
cl 
Q 
| 
ort 
WY Hine 
fom thie eee) 
Ozs0O7 
OWOvee 
HOO 
Be ft 
or) 
oO 
18) 
xi 
-i DY 
an MM 


s: P(D1=1] N4=x) 


MODEL 


NY 
0 
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ODDS 
0.9156 
1.504 


PREDICTED 
FAIL PROBABILITY PROBABILITY 
0.7143 0.7143 
0.8182 0.38182 
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# OF 
mM 
4 


OF 


+ 
SHUecCESS 

19 
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MODEL : 2P(D2=1}D1=x) 


OF # OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS D1 
wOzoU 0.6250 0.511 0 
17 11 0.6071 0.6071 0.435 1 
PARAMETER COEFFICIENT SE COER/SE 
Alpha Oma 11 0.730 0.699 
Beta -0 .076 0.826 -0.091 
Var{ ajJ=0.53 Var{ 8j=0.68 
90% Confidence limits 
LOweL UPDSE 
nouns -0.683 1.705 
Beta -1. 428 1.276 
MODEL : P(D3=1{D2=x) 

* OF # OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS D2 
2D eae ae ap ap as => @» <a» een ee ae eae ee = Seen a 8 o® a a= a ae eeap «a =a a = =» a 

19 4 O27 1483 - 0.7143 0.916 0 
18 4 0.8182 0.8182 1.504 1 
PARAMETER COEFFICIENT Si Beeden 
Alpha 0.916 0.592 1.549 
Beta 0.588 0.810 0.726 
Var{ ,J]=0. 31 Var({ g j=0.52 
90% Confidence limits 
Lower Upper 
nee 0.003 ee 2 
MODEL : P(D4=1]D3=x) 

+ OF # OF OBSERVED PREDICTED LOG 

SUCCESS FAIL PROBABILITY PROBABILITY ODDS D3 
3 5 Os 3750 0.3750 -Q0.511 90 
21 7 0.7500 0.7500 oso tT 
PARAMETER COEFFICIENT SE COEF/SE 
Alpha ~0.511 O3730 ~Qe6ec2 
Beta 1.609 0.851 1.892 
Var{, J=0.53 Var{ 3jJ=0.72 
90% Confidence limits 
Lower ppper 
Alpha =a 10D 0.683 
Beta 02 297 3.001 
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MODEL : P(J1=1|]D4=x) 


# OF # OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS DY 
4 4 0.5000 0.5000 0.000 0 
22 6 0.7857 0.7857 1.299 1 
PARAMETER COEFFICIENT SE CORE ZS. 
Alpha 0.000 0.707 0.000 
Beta 1.299 0.844 1.540 
Var{ ojJ=0.44 Var{ 3 jJ=0.65 
90% Confidence limits 
Lower ota 
meee -1.088 1.088 
Beta -0. 023 2-621 
MODEL ;: P(J2=1|J1=x) 
# OF # OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS J1 
G 6 0.4000 . 0.4000 -0.405 0 
21 5 0.8077 0.8077 1.435 1 
PARAMETER COEFFICIENT Se COEF/SSE 
Alpha -0.405 0.645 -0.628 
Beta 1.880 0.815 Le208 
Var[ , jJ=0.42 Var( g J=0.66 
90% Confidence linits 
Lower Uppes 
eee -1.468 0.658 
Beta 0.508 i ea ar 
MODEL : P(J3=1{[J2=x) 
# OF # OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS J2 
3 8 0.2727 QO22727 -0.981 @) 
ifs] 7 0. 7200 0.7200 0.944 1 
PARAMETER COEFFICIENT SE COEF/SSE 
Beta 1.925 0.810 2.376 
Var{ , J=0.46 Var( g J=0.65 
90% Confidence limits 
LOWer Opper 
Alpha -2.093 0.131 
Beta | 0.603 S247 
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10 DEL : P(J4=1[ J 3=x) 


# OF OBSERVED 


# OF PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS J3 
10 5 0.6667 0.6667 0.693 
17 Q 0. 8095 0.8095 1.447 1 
PARAMETER COEFFICIENT SE COEF/SSE 
Alpha 693 0.548 1.266 
Beta 0.754 0.780 0.966 
Var{ ,]=0.30 Var{ 8 j=0.61 
90% Confidence limits 
Lower 7 eeSs 
pee -0.205 T2591 
Beta -0.527 -035 


MODEL : P(JP=1]{J4=x) 


pe ee Se ee Se ee Se cee 


# OF # OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS J& 
5 3 0.6667 0.6667 0.694 0 
23 4 0.8519 0.8519 1.749 1 
PARAMETER COEFFICIENT SE COBE/SS& 
Alpha 0.694 0.707 0.980 
Beta 1.056 0.891 1.186 
Var{ ,jJ=0.50 Var{ g jJ=0.80 
90% pene eanCe linit 
LO Upper 
aoa -6_ 466 1.854 
Beta -0.411 Ze023 
MODEL : P(F1=1jJF=x) 
OF #% OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS JF 
p 5 0.2857 0.2857 -0.916 G 
25 4 0.8621 0.8621 1.833 1 
PARAMETER COEFFICIENT SE COEF/SSE 
Alpha -0.916 ~837 -1.095 
Beta 2.749 0.995 2.763 
Var{, j=0.70 Var({ gj=1.00 
90% Confidence limits 
Lower Helsis 
a bee -2.288 0.456 
a 1.109 4.389 
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MODEL : P(F2=1}F1=x) 


# OF # OF OBSERVED PREDICTED 1 
SUCCESS PAIL PROBABILITY PROBABILITY ODDS F1 
7 2 0.7778 0.7778 1.050 
20 7 0. 7407 0.7407 1.253 1 
PARAMETER COEFFICIENT SE COBw7s & 
Alpha e253 0.802 1.562 
Beta -0.203 0.914 +-0.222 
Var[ « }]=0.64 Var{ 3 ]=0.84 
90% Confidence limits 
Lower Upper 
Alpha -0. 059 22505 
Beta -1.706 1.300 
MODEL : P(F3=1}F 2=x) 
# OF # OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS eZ 
16 11 0.5926 0.5926 0.375 1 
PARAHKETER COEPFICIENT SE COEBE/SE 
Alpha 0.223 0.671 0.333 
Beta 0.152 0.777 0.195 
Var{ .jJ=0.45 Var( g J=0.60 
90% Confidence limits 
Lower 72535 
Alpha -0. 877 ase S 
Beta -1.118 Nerd 2 
MODEL : P(FM=1{ F 3=x) 
8 OF # OF OBSERVED PREDICTED L 
SUCCESS PAIL PROBABILITY PROBABILITY ODDS P3 
11 uy 0.7333 0.7333 1.012 0 
19 2 0.9048 0.9048 Deon 1 
PARANETER COEFFICIENT SE coErs3° 
Alpha 7.011 0.584 1.733 
Beta 1.241 0.945 1.3171 
Var(, J=0.34 Var({ gJ=0.89 
90% Confidence iimits 
LOWer yPper 
Alpha 0.055 1.967 
Beta -0. 306 2.788 
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MODEL : P(M1=1|FM=x) 


# OF # OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS FM 
5 1 0.3333. 0.8333 1.609 0 
27 3 0.9000 0.9000 2.197 1 
PARAMETER COEFFICIENT SE COEREVSE 
Alpha 1.609 1.095 1.469 
Beta 0.588 e253 0.469 
Varfa J=1.19 Var[ gj=1.57 
90% Confidence linits 
Lower y PSE 
Alpha -0. 180 3.398 
Beta -1. 467 2-643 
MODEL : P(M2=1[M 1=x) 
# OF # OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PR&CBABILITY ODDS M4 
4 0 1.0000 0.9999 2.197 Q 
24 8 0. 7500 6.7500 1.009 1 
PARAMETER COEFFICIENT SE COEF/SSE 
Alpha 2.197 0.000 0.000 
Beta -1.098 0.408 -2.691 
Var{, ]=2.78 Var({ g]=2.94 
90% Confidence limits 
Lower eaogel 
eee -0.537 4.931 
Beta -3.910 1.714 
MODEL =: P(M3=1|M2=x) 
# OF # OF OBSERVED PREDICTED LOG 
SUCCESS FAIL PROBABILITY PROBABILITY ODDS M2 
5 3 0.6250 0.6250 0.511 0 
20 8 0.7742 0.7143 0.916 1 
PARAMETER COEPFICIENT Sis GOES © 
Alpha 0.511 0.730 C.699 
Beta 0.405 0.842 0.482 
Var{ a]=0.53 Var(g J=0.71 
90% Confidence limits 
Lower UPRSE 
Alpha -0. 683 1.705 
Beta -0.977 1.787 
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P(A3=1]A2=x) 
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are dependent on each other at the 10% Significance level 


since the 90% confidence interval for beta does not cover 


zero. 
The following logistic model was also considered. 
meee 
PLY 4) =2IR,} ooo 
l+e 
where Ry is the amount of rainfall in week t and 2S real 1s 


a binary random variable taking the value 1 if week t+l has 
positive rainfall. Tabel 21 shows the alpha and beta coeffi- 
cients, standard errors (ss) and the 90% confidence limits 
for the alpha and beta coefficients. These numbers were 
computed using computer package "BMDPLR". 

The 90% confidence intervals for do not cover zero 
only for the weeks D3 vs D4, JF vs Fl, and M2 vs M3. This 
Suggests that for these weeks the amount of rainfall in the 
previous weeks exceeds the amount of rainfall in the current 
week. Figures 42 and 43 show the plots of estimated alpha 


and beta coefficients by week. 
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EE. Po 


memte 21s RESULT OF LOGISTIC ANALYSIS FOR THE WEEKS 


MODEL =: P(O02=1[{01=x) 


PARAMETER COS FFICIENT SE Soe 
Alpha 0.438 0.409 1.071 
Reta 02027 0.432 0.062 
Var[a ]=0.13 ae 

90% Confidence limit 
Lower gppsz 

Alpha -0. 153 teo2g 

ota -2.771 Zico 
MODEL 3: P (9 3= 1{02= x) 
PARAMETER COEFFICIENT SE eee 
oe -0.735 0.424 =-1.735 
eta 1.067 Ge ler 5 tS 
Varf, J=0. 16 Var{ gj=0.87 
90% Confidence limits 
Lower Upper 
pase -1.391 -0.079 
a -0. 463 2 OST 
MODEL : P(04=1]03=x) 

PARAMETER COESEErCIENT SE COEF/SSE 
Aloha 0.380 0.379 1.002 
Beta 0.317 0.748 0.424 
Var[ aj=0.14 Var{ 8jJ=1.18 

90% SOLE EOS limits 
Lower g°334 
Alpha -0.3384 
Aa -1. 464 098 
MODEL ; P(ON=1]04=x) 

PARAMETER COEFFICIENT SE COEF/SSE 
Alpha -0.192 438 -0.4%39 
Beta 0.636 0.942 0.675 
Var{ ajJ=0.14 Var{ 8j=1. 29 

90% Confidence linits 
LOWer ais 
Alpha -0. 806 Q-822 
Beta =o 227 2.499 


OF 





MODEL : P(N 2=1]N1=x) 
a COEFFICIENT SE 
eee 0.986 0.466 
Beta 0.7483 0.807 
Var{a 1=0.22 Var{ 8]J=0.65 
90% Confidence linits 
Lower Upper 
Alpha 0.21 1.755 
peta -0.579 -065 
MODEL : P(N3=1|N2=x) 
PARAMETER COEFFICIENT SE 
Alpha -0.102 0.422 
Reta 0.284 ae 1 
Varf{a J=0.18 Var{[g j=0.12 
90% Confidence limits 
Lower 2/958 
ee -0.798 O.576 
a -0. 384 852 
MODEL : P{(N4=1] N3=x) 
PARAMETER COEFFICIENT SE 
Aipha 0.595 0.395 
Beta -0.280 0.375 
Varf{a ]=0.16 Var{[s J=0.14 
90% aoe cee linits 
Lower 1384 
ae -0. 061 i 
Beta ~-0. 894 334 
MODEL : P(D1=1| NG=x) 
PARAMETER COEFFICIENT SE 
Alpha 1.106 0.454 
Beta 0.336 0.589 
Vartfa J=0.21 Var{8 J=0.35 
90% Confidence linits 
Lower Fpper 
ae eac 0.354 58 
eta ~0. 634 - 306 
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MODEL 3: P{(D2=1| D1=x) 


=p = ee eee 2 a ae ee ee ee eee ee ee ee Ge 


PARAMETER COEFFICIENT SE COEF/S = 
Alpha 0.294 0.456 0.644 
Beta 0.230 0.450 0.509 
Var{a 3=0.21 Var(8 J=0.20 

90% Confidence Pee 
LOwer i 
Alpha -0.45 P 
Beta -0.59 629685 
MODEL : P(D3=T{ D2=x) 

PARAMETER COEFFICIENT SE COEF/SSE 
Alpha 1.378 0.485 2.838 
Beta -~0.301 0.614  +-0.490 
Var{, ]=0.24 Var(g J=0.38 

90% confidence limits 
Lower ypper 
MCDEL : P(D4= 14 D3 =x) 

PARAMETER COEFFICIENT SE COEF/SSE 
Ipha 0.396 0.482 0.822 
Beta 0.560 0.653 0.858 
Var{a J=0.23 Var{8 J=0.43 

90% Confidence limits 
LOwer queiee 
Alpha -0. 391 1.183 
Beta -0.515 1.635 
MODEL =: P(J1=1] D4U=x) 

PARAMETER COEFFICIENT SE COEF/SSE 
Alpha 0. 132 0.499 0.265 
Beta 1.952 0.988 1.976 
Var{ , ]=0.25 Var{ 2g j=1.07 

90% confidence limits 
Lower aaeSe 
ee -0. 688 0.952 
Beta 0.256 3.648 
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MODEL : P(J2=1]J31=x) 
PARAMETER COEFFICIENT SE 
Alpha 0.299 0.445 
Beta 5232 0.801 
Var{l ajJ=0.20 Var{ 6 J=0.68 
90% Confidence limits 
Lower Gpper 
ele -0. 434 1.032 
Beta -0. 080 2.544 
MODEL : P(J3=1] J2=x) 
PARAMETER COEFFICIENT SE 
Alpha -0.045 0.434 
peta 0.526 0.402 
Var{ ,J=0.19 Var{ 2 J=0. 16 
90% Confidence limits 
LOwer a PpSt 
ae -0. 760 0.670 
Beta -0. 130 1.182 
MODEL =: P(J4=1] J3=x) 
PARAMETER COEFFICIENT SE 
Alpha 0.809 0.472 
Beta 0.394 0.428 
Var{ ojJ=0.22 Var{ §]=0.18 


90% Confadence iiits 


Lower Pees 
Ae 0. O40 1.578 
Beta -0. 302 090 
MODEL ; P(JF=1] J4=x) 
DARAMETER COREPICTENT SE 
Alpha e525 0.567 
Beta -0.137 0.482 
Var{ 4 jJ=0.32 Var[ 6 J=0.23 
90% Confidence limits 
Lower es 
eee US Sem, 2.948 
Beta =Qeo2s 0.650 
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Wovens Dif i=1' JF=x) 
PARAMETER CCl perc Lent SE 
Alpha 0.356 0.505 
Beta 1.162 eee 
Varf{ , J=0.26 Var( gj]=0.4 
90% Confidence limits 
LOwer oss 
Alpha -0. 480 1.192 
Beta 0. 0256 2-298 
MODEL =: P(F2=1] F1=x) 
PARAMETER COEFFICIENT SE 
Alpha 0.986 0.491 
Beta Q.161 0.454 
Var{ a jJ=0.24 Var{ 8 J=0.21 
90% Confidence limits 
LOwer Upver 
eee 0. 183 1.699 
Beta ~0.591 0.913 
MODEL : P(F3=1| F2=x) 
PARAMETER COEFFICIENT SE 
Alpha 0.135 0.444 
Beta O22 0.606 
Var[ , ]=0.20 Var{ g]=0.37 
90% confidence limits 
Lower Koes 
Ae -0. 598 0.868 
Beta -0. 586 1.410 
MODEL : P(FPM=1|] F3=x) 
PARAMETER COEFFICIENT SE 
Alpha 1.4845 0.546 
Beta 0.272 0.578 
Var{ a]=0.30 Var{ 8jJ=0.33 
90% confidence limits 
Lower ees 
Alpha 0. 547 2.343 
Beta -0.670 71.2178 
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MODEL : P(M1=1] FM=x) 
PARAMETER COEFFICIENT SE 
Alpha 1.613 0.634 
Beta 0.637 0.940 
Var{ «j]=0.40 Var[ 8 J=0. 88 
90% Confidence limits 
Lower Upper 
Pope 0.576 2.650 
Beta -0. 701 Zea} 
MODEL : P(M2=1[{M1=x) 
PARAMETER COEFFICIENT SE 
Alpha 1.209 055 25 
Beta 0.083 0.681 
Var{ ,J=0.29 Var{ 2g J=0.46 
90% Confidence linits 
LOwer ae 
nee Ofer 6 2.092 
Beta -1.029 1.795 
MODEL : P(M3=1] 42=x) 
PARAMETER COEFFICIENT SE 
Alpha 0.143 0.467 
Beta : a ae | 0.776 
Var{ aojJ=0.22 Varf 2 j=0.60 
90% Confidence limits 
Lower BPRSS 
Aen -0.626 0.912 
Beta 0.107 2-647 
MODEL : P(M4=1]| 43=x) 
PARAMETER COEFFICIENT SE 
Alpha ng) i | 0.435 
Beta 0.694 
Var({ a jJ=0.19 Var{ 8 J=0.48 
90% Confidence limits 
Lower Upper 
Beta -(). 336 1.936 
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MODEL =: P(A1=1] M4=x) 
PARAMETER Coes he TENT Sz 
Alpha -0.087 0.407 
Beta 0.728 0.466 
Var{ ojJ=0.17 Var{ 8j=0.22 
90% Confidence linits 
LOwer Upper 
wees -0. 763 0.589 
Beta -0. 041 1.497 
MODEL : P(A2=1]A1=x) 
PARAMETER COEFFICIENT SE 
ipha 0.438 0.409 
Beta 0.027 0.432 
Var{ ,j=0.17 Var{ gj=0.19 
90% Confidence limits 
Lower elke 
oleae -0. 238 1.114 
Beta -9. 888 0.742 
MODEL ; P{A3=1| A2=x) 
PARAMETER COEFFICIENT SE 
Alpha -90.735 0.424 
Beta 1.067 0.909 
Var(o jJ=0.18 Var{ 3 j=0. 83 
90% Confidence limits 
LOwer pps 
Seu -1. 431 -0.03 
Beta -0.427 2.561 
MODEL : P(A4&=1|] A3=x) 
PARAMETER COEFFICIENT SE 
Alpha 0.380 0.379 
Beta 0.317 0.748 
Var[ ajJ=0.14 Var{ 8jJ=0.56 
90% Confidence limits 
Lower aeeaa 
Alpha -0.234 0.994 
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MODEL =: DP(AH=1] AY=x) 
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PARAMETER COEFFICIENT SE COEF/SSE 
Alpha -0.192 0.438 -0.439 
Beta 0.636 0.942 0.675 
Var[ oJ=0.19 Var( 8 j=0.89 

90% Confidence limits 
Lower pper 
ae oe . -0.907 0.523 
Beta -Q.9171 2.183 
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APPENDIX C. WEEKLY AMOUNT OF RAINFALL 
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Figure 44. 





YEARS OF THE DATA 


Weekly rainfall in inches for weeks Ol and 02 
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Figure 45. Weekly rainfall in inches for weeks 03 and O4 
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Figure 47. Weekly rainfall in inches for weeks N2 and N3 
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Figure 49. 
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Figure 50. Weekly rainfall in inches for weeks D4 and Jl 


eZ 








WEEKLY AARTNFALL 


6. /\ la Aes \ 
Oi Sr Wi ecoOeme ess SO. 35. 40. 


Cahinor Ome line DATA 


5.7 
et 
F Qu. 
= . 
Li 
- 3. 
oT 
om 
m— 2. 
aa} 

— 
tJ 
~ 12 
= 
6. eesti Get tt Ht T/A te et Se oe Ce ' 3 


QO. Se WOeeeeo wee. 8SU. 35. 40. 


Cat nom Crime WAT A 


Figure 51. Weekly rainfall in inches for weeks J2 and J3 
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Figure 52. Weekly rainfall in inches for weeks J4 and JF 
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Figure 53. Weekly rainfall in inches for weeks Fl and F2 
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Figure 54. Weekly rainfall in inches for weeks F3 and FM 
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Figure 55. Weekly rainfall in inches for weeks Ml and M2 
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Figure 56. WEekly rainfall in inches for weeks M3 and M4 
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Figure 57. Weekly rainfall in inches for weeks Al and A2 
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WEekly rainfall in inches for weeks A3 and A4 
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APPENDIX D 


EXPONENTIAL DISTRIBUTION 


The Exponential Distribution with parameter \ has the 


form 


F(x) = 1 - exp(-Ax) ; Op ex 


= QO ; eu aU) 


It has density function, 


f(x) = dr exp(-rx) 


Suppose X has the exponential density with parameter \. 


Here are some characteristics of X. 


a) The mean ee > ©. O ee i | 
b) The variance : V[X] = We 

std. dev [X] : l1/)d 

uGoef. of variation) - = 1 
c) The Median Xo 5 = 0.693E([X] 
dad) The Lower Quartile: X 45 =e 20 (xX } 
ad) The Upper Quartile: X 46 = 1.386E[(X] 
f) The Skewness : 2 

The Kurtosis : 6 


ALGEBRAIC COMPUTATION OF SKEWNESS AND KURTOSIS FOR EXPONENTIAL 


DESTRIBUTION. 
Exponential distribution with parameter A = 1. 
Density function : £(x) = exp(-x) 


Pi2 





Mean 


Varlance 


A. SKEWNESS 


vol 


Pee KURTOSIS 


Y9 


Eas 1) 


IL 
E((X-1)3] _ a 
tee = B((X-l1)  ) 
J = i (x-1) 2e*ax ; if xXeXax = ki! 
0 0 
= 3! - er + eet - 1x” 
2 


BTGeomy i oe 


J 


oo oO 


If (x-1)%e7*ax; again {| x*e*ax = k:! 


0 0 

4! - 4x7! + 6x7! = 4x7! + 1x91 
9 

= 6 
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SAMPLE PROPERTIES OF SKEWNESS AND KURTOSIS 


Cramer [Ref. 6] gives a discussion of mean and variances 
of the skewness and kurtosis for sampling. In general, the 


mean of Jy and 95 are: 


Elg,] a uy ’ Elg,] = WS 


and the variances are: 


V2 3 2 
Sie p= l2ieawen--24uaul,rou iu 
var(g,] 2° 6 2°35 2a 2. 

4.x nex U5 


+35u5u5+36u> 


Sug Au gl gle usu gugtdug—usugt Ou gush gtlOusus 
ES Ra Tc Mi eM aie) 
a .< Uo 


I 


When the parent population is exponential, 


Elg,! = 27 Elgo] aa 6 


2 





varlg,] = =, varlg,] = 


for the computation. 
In general for the exponential distribution with density 
; -xX TR as 
function f(x) = e “, the kth moment; u, = { Xe “dx = k! 
0 
So, by using this formula and putting it in the equations of 


meariance: 


mi 225 See ey 
var [g,] =) essa Vari[g.,] = 








ESTIMATED AND SIMULATED VALUES FOR SKEWNESS 


AND KURTOSIS FOR SAME SAMPLE SIZES. 


Table 22 
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# OF YEARS 
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